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ABSTRACT 

The inhomogeneous cooling flow scenario predicts the existence of large quantities of gas in massive 
elliptical galaxies, groups, and clusters that have cooled and dropped out of the flow. Using spatially 
resolved, deprojected X-ray spectra from the ROSAT PSPC we have detected strong absorption over 
energies 0.4 - 0.8 keV intrinsic to the central ^ 1' of the galaxy, NGC 1399, the group, NGC 5044, 
and the cluster, A1795. These systems have amongst the largest nearby cooling flows in their respective 
classes and low Galactic columns. Since no excess absorption is indicated for energies below ~ 0.4 keV 
the most reasonable model for the absorber is warm, coUisionally ionized gas with T — 10^~® K where 
ionized states of oxygen provide most of the absorption. Attributing the absorption only to ionized gas 
reconciles the large columns of cold H and He inferred from Einstein and ASCA with the lack of such 
columns inferred from ROSAT, and also is consistent with the negligible atomic and molecular H inferred 
from Hi and CO observations of cooling flows. The prediction of warm ionized gas as the product of 
mass drop-out in these and other cooling flows can be verified by Chandra, XMM, and ASTRO-E. 

Subject headings: cooling flows - intergalactic medium - X-rays: galaxies 



1. INTRODUCTION 

The inhomogeneous cooling flow scenario (e.g., Fabian 
1994) is often invoked to interpret the X-ray observations 
of massive elliptical galaxies, groups, and clusters. The 
key prediction of this scenario is the existence of large 
quantities of gas that have cooled out of the hot phase 
and dropped out of the flow. The only evidence for large 
amounts of mass drop-out arises from the excess soft X-ray 
absorption from cold gas found for many cooling flows es- 
pecially from spectral analysis of Einstein and ASCA data 
(e.g.. White et al 1991; Fabian et al 1994; Buote & Fabian 
1998; Buote 1999, 2000a). 

This interpretation is highly controversial because for 
systems with low Galactic columns no excess absorption 
from cold gas is ever found with the ROSAT PSPC which 
should be more sensitive because of its softer bandpass, 
0.1-2.4 keV (e.g., David et al 1994; Jones et al 1997; Briel 
& Henry 1996). Furthermore, the large intrinsic columns 
of cold H indicated by the Einstein and ASCA data are in 
embarrassing disagreement with H i and CO observations 
(e.g., Bregman, Hogg, & Roberts 1992; O'Dea et al 1994). 

We have found new evidence for absorption in the PSPC 
data during our investigation of the radial metallicity pro- 
files of several of the brightest, nearby cooling flows (Buote 
2000b, hereafter PAPER2). In this Letter we present ab- 
sorption proflles for a subset of these cooling flows obtained 
from analysis of the deprojected PSPC spectra and briefly 
compare these results to single-aperture ASCA spectra. 
We focus on the galaxy, NGC 1399, the group, NGC 5044, 
and the cluster, A1795 because these systems have (1) 
low Galactic iVn which facilitates analysis of any intrinsic 
absorption, (2) the most signiflcant excess columns mea- 
sured from two-component spectral models with ASCA, 
and (3) the most clearly signiflcant differences in A^h ob- 
tained from PSPC data when energies below ~ 0.5 keV are 



included/excluded from analysis. The temperature and 
metallicity profiles, as well as results for a few other sys- 
tems not satisfying all the above conditions, will be given 
in PAPER2 and Buote (2000c, hereafter PAPER3). 

2. SPATIALLY RESOLVED ROSAT SPECTRA 

We obtained archival ROSAT PSPC observations of 
NGC 1399, NGC 5044, and A1795 and reduced the data 
as described in PAPER2. Here we mention that particu- 
lar attention was given to excluding any significant fluc- 
tuations in the light curves since contamination from so- 
lar emission would preferentially appear in the low-energy 
channels. We also developed our own software to properly 
scale the background spectra to the source positions. 

For each object we extracted spectra in concentric cir- 
cular annuli located at the X-ray centroid such that for 
each annulus the width was > 1' and the background- 
subtracted counts was larger than some value chosen to 
minimize uncertainties on the spectral parameters for each 
system while maintaining as many annuli as possible. Data 
with energies < 0.2 keV were excluded to insure that the 
PSF was < 1' FWHM. For our on-axis sources ~ 99% of 
the PSF at 0.2 keV is contained within R = 1'. 

We deproject the data following McLaughhn (1999) who 
elaborates on the original paper by Fabian et al (1981). 
The spectral fitting is performed with XSPEC vlO.O using 
the MEKAL code to model the (single-phase) hot plasma 
emission. We used the photoelectric absorption cross sec- 
tions of Baluciiiska-Church & McCammon (1992). Al- 
though Arabadjis & Bregman (1999) point out that the 
He cross section at 0.15 keV is in error by 13%, since we 
analyze i? > 0.2 keV we find that our fits do not change 
when using the Morrison & McCammon (1983) cross sec- 
tions which have the correct He value. Further details of 
our deprojection analysis are discussed in PAPER2. 
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2.1. Intrinsic Absorption from Cold Gas 

No excess absorption from cold gas is found in cooling 
flows with the PSPC when energies down to ~ 0.1-0.2 
keV are included, whereas significant absorption (similar 
to that obtained by ASCA) is found when the PSPC spec- 
trum is restricted to energies above ^ 0.5 keV (Allen & 
Fabian 1997; Buote 1999). Presently, the only viable ex- 
planation of why large excess absorbing columns are not 
inferred when energies down to ~ 0.1-0.2 keV are included 
is that the standard foreground screen model systemati- 
cally underestimates the true column intrinsic to the cool- 
ing flow (Allen & Fabian 1997; Sarazin, Wise, & Marke- 
vitch 1998). Allen & Fabian show in their Figure 9 that 
a simulated PSPC cluster spectrum modified by both a 
foreground column of iVn = 10^° cm~^ plus an intrinsic 
column of TVh = 10^^ cm~^ with covering fraction 0.5 will 
always yield a value similar to the foreground column if fit- 
ted only with a standard foreground model. However, such 
two-component models usually represent the projection of 
extended gas (with no intrinsic absorption) with centrally 
concentrated, intrinsically absorbed gas, and thus the suit- 
ability of these models can be tested via deprojection. 

NGC 1399: R = 0-1 arcmin 
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Fig. 1. — ROSAT PSPC spectrum and projected model (see 
text). Statistical errors on the model are similar to those shown for 
the data. 

In Figure |l| we display the PSPC spectrum for NGC 
1399 within R = 1' (2D) and show the model representing 
the projection of gas from r > 1' (3D) into the R — V 
aperture modified by a foreground Galactic absorber with 
solar abundances. The deprojection separates the fore- 
ground component (i.e. the projected model in Fig [l]) 
from the intrinsic gas and shows that the intrinsic com- 
ponent cannot be absorbed very differently from the outer 
gas since their spectral shapes below ^0.5 keV are so sim- 
ilar; e.g., for an intrinsic absorber with covering factor of 
0.5 we obtain a best-fit AiVn = and AiVn < A^h'"' at 
> 90% confidence. Since we obtain analogous results for 
NGC 5044 and A1795, we conclude that the deprojected 
PSPC spectra for each system rule out absorption by large 
amounts of intrinsic cold gas (especially H and He) as a vi- 
able explanation of the different excess absorbing columns 
obtained from analyses in different X-ray bandpasses. 

2.2. Oxygen Edge 

A solution to this problem becomes apparent upon ex- 
amination of how Nn varies with the lower energy limit in 
the centers of cooling flows. For example, in the central 1' 



of NGC 1399 we find that Ar ~ N^""^ when E^i^ 0.2- 
0.3 keV, though when i^min ^ 0.4 keV A'h increases to 
~ 2N^'^^. A dramatic change occurs for i?niin ^ 0.5 keV 
where A^h increases to several times A^^^' , although for 
larger i^min A^h remains nearly constant. In Figure ^ we 
plot A^H {R) for E^in = 0.2, 0.5 keV for NGC 1399, 5044, 
and A1795. The profiles radically differ for each Eynin- 
(R) - A^ciai ^^^^ ^^.^ ^ 0.2 keV, but A^h (R) in- 
creases significantly at small R when Emin = 0.5 keV such 
that A^H > Ag^' at > 90% confidence within R = V . (The 
uncertainties are computed from 100 Monte Carlo realiza- 
tions of the best-fitting models in each 2D annulus. Our 
simple procedure for defining confidence limits (see PA- 
PER2) occasionally results in the best-fitting value lying 
outside the la errors.) 

Since Ah (i?min) ~ const for i^min > 0.5 keV the por- 
tion of the spectrum responsible for the absorption must 
be 0.4-0.7 keV considering the PSPC resolution. The 
most important absorption (and emission) features over 
these energies are due to oxygen. In Figure || we also 
show the deprojection results for i^min — 0.2 keV includ- 
ing an edge at 0.532 keV (rest frame) corresponding to 
cold oxygen (Oi). In this case A'h (R) is very similar to 
the previous fits with Smin — 0.2 keV and no edge (i.e., 
similar to Galactic), and the edge optical depth profiles, 
t{R), closely resemble A^h {R) when E'min = 0.5 keV. 

The large optical depths at small radii and the improve- 
ments in the fits when adding the edge are highly signifi- 
cant. In the R = 1' bins the 95% confidence lower limits 
on T are 0.97, 0.59, and 0.21 respectively for NGC 1399, 
5044, and A1795. Furthermore, 99 out of 100 Monte Carlo 
simulations predict r larger than 0.31, 0.10, and 0.17 re- 
spectively; i.e., T > is highly significant for these systems 
in the inner radial bin. The fits are improved substantially 
when the edge is introduced in the central bin of each sys- 
tem - see Table |l| where we list the values of x^, degrees of 
freedom (dof), and null hypothesis probability (P). (Note 
that for NGC 1399 the edge improves the fits significantly 
more than by letting the oxygen abundance in the hot gas 
go to zero: i.e., improves only to 95.9 as opposed to 
83.7 for the edge.) We emphasize that for the large values 
of r ^ 1 obtained in the central regions the absorption 
from the edge affects a large energy range comparable to 
the resolution of the PSPC over ~ 0.5-0.7 keV; e.g., the 
Oi edge absorbs 25% of the flux at 0.8 keV for r = 1. 

For most radii the constraints on the edge energy are 
not very precise which is why we fixed the edge energy in 
our analysis. If we allow the edge energy to be a free pa- 
rameter in the central radial bin we obtain 0.51+5]:[J^ keV, 
0.511° °^ keV, and 0.561° °^ keV (90% confidence) for the 
edge energies of NGC 1399, 5044, and A1795. These con- 
straints are consistent with the lower ionization states of 
oxygen but not edges from the highest states Ovi-viii. Due 
to the limited resolution we can add additional edges to 
share the r obtained for the O i edges, although even with 
a two-edge model a significant r cannot be obtained for 
edge energies above 0.65 keV corresponding to ^^Ovi. 

(We mention that the models for NGC 5044 with i^min = 
0.2 keV not only imply A'h < A^*^' for small R but also 
predict very large metallicities that are inconsistent with 
the ASCA data. If A'h is fixed to N^'^^ then consistent 
metallicities are obtained as are t{R) and improve- 
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Fig. 2. — (Upper panels) Radial column density profiles and Icr errors obtained from deprojection analysis of PSPC spectra for fits over 
0.2-2.2 keV (solid line), 0.5-2.2 keV (dashed), and over 0.2-2.2 keV but including an O edge (dotted). (Lower panels) Optical depth of 
0.532 keV O edge (at redshift of object); dashed diamonds for NGC 5044 represent case where Nn = A^^^' . In the outermost bins we set 
A^H = N§^^ for 

^miii — 0.5 keV and t = because of weak constraints in some cases. Radial units in kpc on top axis assume distances of 
18, 38, and 240 Mpc respectively for NGC 1399, 5044, and A1795. 



Table 1 

X^/dof/P for Model Fits to Deprojected PSPC Spectra within r = 1' 



NGC 1399 NGC 5044 A1795 

Free A^h A^h = A^g"-' 

no edge 104.9/96/0.25 138.6/117/8. 5e- 2 191.2/118/2.3e-5 253.2/182/3. 8e-4 
w/ edge 83.7/95/0.79 126.4/116/0.24 156.5/117/8.7e-3 212. 0/181/5. 7e-2 



ments very similar to the free A'h case - see Figure g and 
Table |l|; see also PAPER2 and PAPERS.) 

3. SINGLE-APERTURE ASCA SPECTRA 

The intrinsic oxygen absorption indicated by the PSPC 
data is most significant in the central 1' which is much 
smaller than the width of the ASCA PSF. In addition, 
since the ASCA SIS is Hmited to > 0.5 keV, and the ef- 
ficiency near 0.5 keV is limited due to instrumental oxygen 
absorption, it cannot be expected that ASCA can distin- 
guish an oxygen edge from an absorber with solar abun- 
dances. Nevertheless, we briefly summarize the results ob- 
tained when adding an oxygen edge to the ASCA data. 

Previously we (Buote 1999) have fitted two-temperature 
models to the accumulated ASCA SIS and GIS data 



within i? « 5' of NGC 1399 and 5044 and obtained 
= 49;g, 25^6 X 10^° cm~2 respectively for the cooler 
temperature components. Examination of Figure ^ shows 
that these columns are ^ 2 times the values in the cen- 
tral arcminutc and are consistent with the total columns 
within i? « 5' obtained from the PSPC for i?min = 0.5 
keV. If instead we add an Oi edge to the cooler ASCA 
model components (while keeping A^h fixed to Galactic 
on each component) we obtain (1) fits of comparable qual- 
ity (slightly better) to the original models and (2) opti- 
cal depths that agree with the cumulative values obtained 
from the PSPC. 

For A1795 we have re-analyzed the ASCA data within 
R Ki A' and obtained results very similar to Fabian et al 
(1994) for two-temperature models. As above, we obtain 
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fits with an oxygen edge supplying the excess absorption 
that are as good as those obtained with an excess absorber 
having solar abundances. Unlike the ellipticals and groups, 
for the cluster A1795 the cooler temperature component 
contributes only ~ 1/7 to the emission measure, and thus 
its fitted parameters are not nearly as well constrained 
as for NGC 1399 and 5044; i.e., although the best-fitting 
T = 2.5, the 90% lower limit is 0.6 which is comparable to 
the total optical depth inferred from the PSPC (Fig ^). 

Hence, the oxygen edge provides as good or better de- 
scription of the excess absorption inferred from multitem- 
perature models of ASCA data within the central few ar- 
cminutes of NGC 1399, 5044, and A1795 as a cold absorber 
with solar abundances, and yields optical depths that are 
consistent with those obtained with the PSPC data. 

4. DISCUSSION 

Using spatially resolved, deprojected X-ray spectra from 
the ROSAT PSPC we have detected strong oxygen absorp- 
tion intrinsic to the central ^ 1' of the elliptical galaxy, 
NGC 1399, the group, NGC 5044, and the cluster, A1795 
which are amongst the largest nearby cooling flows in their 
respective classes and have low Galactic columns (e.g., 
Fabian et al 1994; Buote & Fabian 1998; Buote 1999, 
2000a). Modeling the oxygen absorption with an edge 
(rest frame E = 0.532 keV) produces the necessary ab- 
sorption in both the PSPC and ASCA data for E ^ 0.5 
keV without violating the PSPC constraints over 0.2 ^ 0.3 
keV for which no significant excess absorption is indicated. 
This reconciles many reported discrepancies between ab- 
sorbing columns inferred from ROSAT with those obtained 
from ASCA and other instruments with bandpasses above 
~ 0.5 keV. Moreover, there is no need for large absorbing 
columns of cold H which are known to be very inconsistent 
with the negligible atomic and molecular H measured in 
cooling flows (e.g., Bregman et al 1992; O'Dea et al 1994). 

Since no excess absorption is detected over 0.2 ^ 0.3 
keV, and since the edge energies allow for ionized states 
of oxygen, the most reasonable model for the absorber 
is warm, coUisionally ionized gas with T = 10^~^ K. At 
these temperatures the majority of the absorption arises 
from oxygen but with a sizeable contribution from ionized 
carbon and nitrogen (see PAPER3). 

Arnaud & Mushotzky (1998) have reported the detec- 
tion of an O I edge in the Perseus cluster using BBXRT 



data. Their cooling flow model (which has stronger oxy- 
gen lines than two-temperature models) requires an edge 
energy that is consistent only with cold oxygen, and thus 
warm ionized gas could not be responsible for the absorp- 
tion. However, it was not shown that multiple edges from 
different ionization states are inconsistent with the data, 
and thus we expect that a multi-edge model would fit as 
well as a single Oi edge since the BBXRT has lower en- 
ergy resolution than ASCA. (We do not find any evidence 
for excess oxygen absorption in Perseus from the PSPC 
data probably because of the large Galactic column, but 
T ^ 2 for the O i edge obtained by Arnaud & Mushotzky 
is similar to what we found for A1795 with ASCA.) 

The warm ionized gas implied for NGC 1399, 5044, and 
A1795 (and possibly Perseus) could be the much sought- 
after mass that has dropped out of the cooling flows. In 
PAPER3 we show that the absorber masses indicated by 
the oxygen edge optical depths are consistent with the 
mass expected to have been deposited by the cooling flow 
over the lifetime of the system, while the emission ex- 
pected from the warm gas does not violate published op- 
tical and UV constraints. Understanding the details of 
how warm ionized gas is maintained and supported in such 
large quantities is a serious theoretical challenge. 

Since detection of the warm ionized gas at other wave- 
lengths is difficult, e.g., background QSOs have not been 
found within i? = 1' of strong cooling fiows (Miller, Breg- 
man, & Knezek 2000), confirmation will probably have to 
be realized in the X-ray band. Fortunately, the CCDs of 
the Chandra and XMM missions have the combined spa- 
tial and spectral resolution and include energies down to 

0.1 keV to verify our prediction of warm ionized gas in 
these and other cooling fiows. The ^ 10 eV resolution of 
the ASTRO-E XRS will place strong constraints on the al- 
lowed ionization states of oxygen. Hence, data from these 
new missions will elucidate the properties of the warm ion- 
ized gas and thus the role of cooling flows in the formation 
and evolution of ellipticals, groups, and clusters. 

I thank M. Bolte and the referee for comments on the 
manuscript. Support for this work was provided by NASA 
through Chandra Fellowship grant PF8- 10001 awarded by 
the Chandra Science Center, which is operated by the 
Smithsonian Astrophysical Observatory for NASA under 
contract NAS8-39073. 
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